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The interior of a vacuum bubble in de Sitter space may give an open universe with sufficient 
homogeneity to agree with observations. Here, previous work by Bucher, Goldhaber and Turok 
is extended to describe a thin bubble wall with nonzero radius and energy difference across the 
wall. The vacuum modes present before formation of the bubble propagate into the interior of the 
open universe and the power spectrum of the resulting gauge invariant gravitational potential is 
calculated. It appears to become scale invariant on small scales, with onset at about the same scale 
as that found in the zero radius case. There is sensitivity to the radius and energy difference at 
large scales, but it is expected that they cannot be strongly constrained because of cosmic variance. 
The prediction of a scale invariant spectrum seems to be robust with respect to variation of these 
parameters at small scales, and apparently is a generic feature of the contribution of these modes 
for these thin wall models. 
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I. INTRODUCTION AND BACKGROUND 



Inflationary models usually predict = 1 but it has recently been shown they also allow for the possibility 
of open and closed § universes. The current models of open universes are either based on nucleation of 

vacuum bubbles^,^ or on open universes which are also compact Q . This paper extends work on models based on 
bubble nucleation. 

Bubble nucleation occurs when a field configuration trapped in a false vacuum tunnels through a barrier to the 
true vacuum. The tunneling only happens over a finite region of space, so that a bubble of true vacuum forms, 
with a bubble wall (a domain wall) interpolating between the false vacuum outside and the true vacuum inside. The 
thickness of the bubble wall and the size of the bubble are determined by parameters in the theory. If there is more 
than one field, the bubble wall interpolates between true and false vacuum values of all the fields in the theory. 

The interior of a vacuum bubble is an open universe [||J|] . If there are random initial conditions before nucleation, 
the open universe in the bubble interior has primordial inhomogeneities, which may be erased by a subsequent phase 
of inflation. However, the amount of inflation required to obtain sufficient homogeneity to agree with observation 
will also make the universe approximately flat, with a critical density f2 « 1. Recently proposed |^-|^ bubble models 
accomplish homogeneity without flatness by having three stages. First there is a phase of ordinary inflation, driving 

— > 1 and erasing initial inhomogeneities. Then a bubble is nucleated, resetting 17 ~ in the interior of the bubble. 
A period of slow roll inflation then follows, solely for the purpose of driving fl up to some desired value slightly less 
than one. At the end of this second stage of inflation, the usual inflationary scenario for reheating takes over. The 
open universe is the interior of the forward light cone of the center of the bubble. The bubble wall, exterior to this 
light cone by definition^, has moved out to infinity before some time T in the past. 

Modifications using two fields [§j5[0] can result in universes with variable values of One field provides infiation 
and the other field resets 17 ~ through bubble nucleation at some point. The advantage of this is that more natural 
(e.g. polynomial) potentials can be used for the fields. (The one field models need more complicated potentials 
to resemble the calculable toy models, although these might occur in a supergravity model.) Some care must be 
taken to avoid large inhomogeneities due to the change of natural time coordinate after tunneling (the synchronicity 
problem). In addition, predictability can be limited because the amount of infiation, depending on one of the fields, 
may not necessarily be correlated with the random tunneling time of the other field, producing a range of different 
f2's. However, there are modifications which can give more predictability As mentioned in ||^, much of the one 

field model analysis carries over to two field models. 

The three stage scenario has some general properties [1): Before the bubble nucleates, the first stage of infiation 
drives spacetime to the de Sitter invariant vacuum, the Bunch-Davies vacuum. The Bunch-Davies vacuum was argued 
in 1^ to be the natural vacuum because it has the right symmetries and it goes over to the Minkowski vacuum at 
short distances. This provides some robustness with respect to initial conditions. The recurrence of inflation after 
tunneling means the bubble must tunnel from de Sitter space to de Sitter space. The slow roll in the second phase of 
inflation means that the effective mass after tunneling should be approximately ~ 0. The second stage of inflation 
must be fine tuned (i.e. the parameters in the potential must be fine timed) within a few percent in order to get 
17 - .1 - .9. 

It is of interest to know how the particulars of the bubble formation affect the resulting density perturbations. 
Unlike ordinary inflation, here some fine tuning of the parameters in the potential is required to fix 17, so it is possible 
that other properties of the resulting theory are sensitive to variation of potential parameters. The focus here is on 
the imprint of the bubble wall on the vacuum density perturbations. There are also perturbations arising from the 
fluctuations of the bubble wall itself P,[7|,p^|-p^ which must be included. 

In order to calculate the density perturbations one does the following: The density perturbations start out (since 
an initial long period of inflation is assumed) as Bunch-Davies vacuum modes of the false vacuum corresponding to a 
field of some mass M . Tunneling is formally a process in Euclidean time, so these modes are analytically continued 
to Euclidean time to give initial conditions. The bubble wall is a boundary between two de Sitter spaces, which 
will be called the false and true vacua from here on (the final space is not actually a vacuum but the slow rolling 
second stage of inflation, it will be called the true vacuum here since it is the endpoint of the tunneling). The initial 
false vacuum fluctuations propagate through the bubble wall and are matched onto a configuration of true vacuum 
modes using Bogoliubov coefficients. At the nucleation time, the whole system is analytically continued to Lorentzian 
signature. Using the natural coordinate system for matching across the wall means that the resulting wavefuction 
must be normalized by hand at some point and can be done here. The perturbations, once inside the forward light 



* although for a thick wall there may be some overlap 
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cone of the bubble (which requires another analytic continuation, since the natural coordinate system is singular on 
the light cone), are related to perturb ations of the gauge invariant gravitational potential, $. The power spectrum of 
eq. ( |5.9| ) along with eqns.(|46,4.7), is the main result of this paper. It can be combined with contributions from 
the wall fluctuations and some model of matter content in the universe (e.g., cold dark matter) to get temperature 
fluctuations to compare to the CMB. 

The formalism for fluctuations around a tunneling background has been developed in many papers, see and 
references therein. Bogoliubov coefRcients for specific thin and thick bubble walls (both designed to be tractable and 
yet have enough parameters to be able to test for robustness), were calculated in Minkowski space in In these 

the fluctuations were massive in both the true and the false vacuum. 

The inclusion of gravity and the true vacuum mass m? ~ add extra complications. In Q], the wall was taken to 
be a step function, and the radius of the bubble at nucleation, R, was taken to be zero. The power spectrum was 
calculated for a field with initial false vacuum mass = 2H^, and the metric perturbations corresponding to the 
field fluctuations were identified. They also found the matching conditions for a masslcss field across the light cone, 
which is subtle due to only half of the modes of a massless field propagating into the light cone. A second paper, [ pO[ , 
generalized this to a field with arbitrary mass in the false vacuum, but still with radius zero. In |18|, the quantum 
state was found for a thin wall bubble with gravity when the fluctuations only have mass on the wall itself. 

Effects of perturbations of the wall itself have also been studied in |0,|io|,^ , building on work on bubble fluctuations 
by The temperature fluctuations predicted in the CMB from Bunch-Davies vacuum fluctuations (without a 

bubble) have been calculated in |^ and compared to the CMB measurements and predictions of the conformal ]2^ ] 
vacuum. 

In this paper, the thin wall approximation is retained, but the bubble wall is given a finite radius R, for an initial false 
vacuum mass conformal or larger, > 2H^. Some of the other assumptions used previously are also carried over: 
that the era preceeding tunneling has resulted in a homogeneous Bunch-Davies vacuum as initial conditions and that 
bubble nucleation is rare enough that no other bubbles are nearby. First the background solutions, approximations and 
fluctuation modes are reviewed. After finding Bogoliubov coefficients, the wavefunction is normalized and continued 
into the forward light cone of the bubble interior (the open universe) . Then the power spectrum for the gauge invariant 
gravitational potential is calculated for this specific background using the results of |2j and compared to earlier cases. 
These results must be combined with those for fluctuations of the wall itself and some model for structure formation 
{e.g., cold dark matter) in order to make a prediction for experiments. 

Note Added: While this paper was being prepared for submission, the preprint |24[ appeared, which independently 
considers some of the same issues. Some puzzles pointed out in this paper are addressed therein; the final version of 
this paper has been updated to refer to these where appropriate. 



II. BACKGROUND SOLUTION, COORDINATE SYSTEMS, AND APPROXIMATIONS 

The background bubble solution for the fluctuations is reviewed first The bubble is described to first 

approximation by classical evolution in Euclidean space. After nucleation, the system is rotated to Lorentzian signature 
and continues evolving classically . 

The action for a single scalar field plus gravity is 

S = [d^x V^ilg^^'d^^d^cj, - V{4,) - {IGTTGy^R] (2.1) 
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with metric g^^,, curvature scalar R, and potential V{4>). In the flat space case, a solution respecting S0(4) symmetry 
has minimum Euclidean action, so when gravity is included, one looks for a solution depending only upon the 
rotationally invariant radial coordinate a. This has not been proven to be a property of the minimum action solution 
for the case with gravity, but none with smaller action are known. 

With gravity, the Euclideanized metric for tunneling is written in terms of , 

^ dCT^ + 6^(cr)[dT| +C0sV£;rff7(2)] , (2.2) 

with dfl(2) = de^ + sin^ Odcj)'^ . 

The Euclidean equations of motion are then 

and 
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b'{a)Y _ 1 SttG 



(2.4) 



The boundary conditions are 



b{a) J P{a) i 

b{a = 0) = b{(7 = amax) = , , . 

</,'(a = O) = 0'(a = a™a,)=O. ^^'^^ 



where the turning point in the tunneUng solution, amax, is infinite for the flat space case, but becomes finite once 
gravity is included p9| , p6[ . For « constant, define = 87rGV{(j))/3, so that for de Sitter space, b{a) = shi[Ha\. 
For fiat space, b{a) = a. 

For many of the calculations it will be useful to switch to a second coordinate system, defined by tanhii = cos[H(j\. 
In terms of these coordinates, the Euclidean line element is: 

ds^ = a'^{u){du'^ + dT% + cos^ TEdn(^2)) , (2.6) 

for de Sitter space, 

H ^ cosh u 

As one approaches the light cone, cr — s- 0, and u — > oo 



The equations of motion (2.3 2.4) have not been solved in general. A few different approximations are possible, 
whose validity depends on parameters in 1^(0). A standard approximation is the thin wall approximation, where the 
energy difference between the true and false vacuum, — Hp — H^, is very small ||^. The bubble wall, the boundary 
between the two vacua, is then located approximately at one value of b{a) = R. In order to consistently neglect r^, 
the width of the wall, it must be small enough compared to other parameters in the theory, 

rb < i? , rfc < , rb < (i7| + H^)-^^^ = A'^ . (2.8) 

In order to verify these conditions, the radius R in the presence of gravity can be calculated using the radius R in the 
absence of gravity, A, and e ||28|| . 

Another possible approximation is to neglect gravity. The validity of this is measured by the second term in equation 



(2^), since for flat Minkowski space b'{cr) = 1. 

In j^] a 'new thin wall' approximation was made. They assumed that R < H^^, the wall is thin, the potential is 
slowly varying after nucleation (Mtme ~ 0), and that the bubble nucleates at about the same energy as that of the 
false vacuum (e « 0). This allows one to treat the background as fixed and to neglect gravity, which they verified is 
consistent. The energies of the initial and final states are approximately degenerate. Hp « Ht, and i? « 0. After 
tunneling, the potential is flat enough for slow roll inflation, so that it is flat enough for the field (f> to be considered 
constant. Thus expressions for tunneling from de Sitter space to de Sitter space can be used. 

The background for fluctuations in this paper is the thin wall solution found by Eq] , which, in terms of the coordinate 



u above, is 



= i ^ u>UR (2-9) 

Hp [cosh( u—ur) cosh UR-\-sinh.{u-~-UR)'\/ cosh^ ur — {Ht / Hp)'^^ 



Mu)^{t l^'V (2-10) 

This solution requires ur > 0. The interior a(u) can be recognized as {Ht cosh{u + S))^^, the metric for a shifted 
value of u, so that the metrics are matched at the same value of a{u), which is the radius 

1 1 

z = I (2 11) 

Hp cosh UR Ht cosh(-uij: + 6) 



The coordinate shift 6 obeys 



(1 - VI - {HpRmi + vi - (htry 



[HtR){HpR) 



(2.12) 
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Setting the Ht = Hp takes 5 — > 0. Taking the matching radius i? means (from eq. ( pi) )) Ht ^ Hp and 
ur oo. 

Once the bubble is nucleated, one rotates to a Lorentzian signature to continue the evolution. Usually the nucleation 
time is chosen at = 0. The coordinate system for the Euclidean description of the bubble can be analytically 
continued, 



ds^ = da^ + b^{a)[-dT^ + cosh'' Td^l(2)] 
= Wrtsh^idu^ - + cosh^ rdl7(2)] 



(2.13) 



The equations of motion for the a dependent background remain unchanged. 

On the light cone of the bubble center, b{a) = (u = oo), so there is a coordinate singularity. One can continue to 
the interior of the light cone by rotating, a = ~iT, t = x ^ i7r/2, b — —ic, to the coordinate system with line element 



ds^ = -dT^ + c\T)[dx^ + smhxdn^2)] ■ (2.14) 



The equations of motion are now 



c(T) + 3fi0(r) = -f 



(2.15) 



For de Sitter space, the background solution except at the wall (which is not inside the interior of the light cone for 
a thin wall) is c{T) = H-^ smh[HT]. 

Classically, then, the following picture arises there are two regions of space considered here, outside the light 
cone of the bubble center (which will be called region II, as was done in |^) and inside the light cone of the bubble 
center (region !).[] For the thin and new thin wall approximation, the bubble wall is approximately located at just 
one value of b{a) = R. The wall moves exterior to the light cone of the bubble's center, approaching it on the curve 
6(cr)^ = R^. (For R = the bubble wall is the light cone.) The bubble wall has no intersection with any of the curves 
in region I. There may be some overlap with region I for a wall that isn't thin. The energy difference between the 
true and false vacuum goes into accelerating the wall (which as it increases in size, due to surface tension, contains 
more energy). 

In region I, inside the future light cone, the S0(3,l) invariant line element behaves just like time in the equations 
of motion. The field (j) is constant on lines of constant T and thus region I can be viewed as an open Freedman- 
Robertson Walker universe with scale factor c(T). The constant T surfaces are not Cauchy surfaces, however, which 
will be important later on when the spectrum is considered. Since a thin bubble wall is confined to region II, it occurs 
"before" T = for the open universe. 

There are consistency requirements for parameters in a given background theory, some of which are automatically 
satisfied in the thin wall approximation. The action of the bubble should be large in order for the semiclassical 
approximation to be valid, and can be calculated directly from the action in the absence of gravity p8| . A large 
action makes the nucleation rate small, another requirement for the scenario here, so that tunneling occurs only after 
inflation has driven the initial conditions into the Bunch-Davies vacuum. The surface tension controls how much the 
bubble wall remains close to its classical configuration, and thus it must be large for the classical solution to be a good 
saddle point. Increasing the surface tension increases the action of the bubble, so these two conditions are related. 
If the false vacuum mass is too small, the lowest action solution is the Hawking Moss solution, ||2^ corresponding to 
stochastic fluctuations over the top of the barrier rather than tunneling through it. This gives unacceptable density 
perturbations, ||§,|3,|o|. 

It is not clear how difficult it is to construct potentials which can give rise to the general thin wall background 
used here. Simple (j)'^ potentials in one field models tend to result in the Hawking Moss instanton when the slow 
roll condition in the true vacuum is combined with thin wall conditions. A linear potential in both the true and false 
vacuum (triangular barrier) apparently forces the energy difference to be negligible. In addition to finding out which 
parameter values are consistent with the thin wall approximation, it is particularly interesting to see what potentials 
arise naturally in particle physics models. 



^To cover de Sitter space requires another region antipodal to I, but it will only be relevant here when verifying the wavefunction 
normalization inside region I, summarized in section 5. 
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III. FLUCTUATIONS AROUND THE BACKGROUND SOLUTION 



With the background solution of the previous section, eq. ( p.9D , one can now solve for basis functions in both the 
true and false vacuum and match the false vacuum fluctuations across the wall. (A rigorous foundation for this using 
the WKB formalism is found in See also [^-16, ij] and references therein.) 

The potential right before nucleating the bubble is assumed to be a deep well, and after nucleating practically flat, 
for the second bout of inflation. Schematically, 

For a thin wall, the matching is in region II. The discussion of basis functions as in |^] is reviewed first. The region 
II fluctuations around the classical solution will be given for Lorentzian signature and then rotated to Euclidean 
signature. They obey (since the background is constant except at the wall) 

{O-^-^J^U=4>JSc^^{O~fn'{<7m^0 (3.2) 

The mass only depends upon a because the background solution {(/);, (ct), b{a)} only depends upon a. So with a step 
bmction at b{a) = R, the boundary conditions are 

2. . ( ]VP , R< b{a) <oo ,„ 
('^^ = I , < b{a) < R (3-3) 



Expanding cq. (3^) in terms of (a, r, 6^ <p)^ and the background de Sitter metric, 

[ +3H cot Ha- _^^( +2tanhT— + — ^) - m2(a)]0(a, r, 0, 0) = 0. (3.4) 

oa^ OCT sm Ha ot^ ot cosh r 

The eigenfunctions of the Laplacian on the 3d hyperboloid satisfy 

g2 Q ^2 

(— + 2tanhr— + -^)Yphn{T,0,4>) = -{1 + p^)Ypi„^{T,e,^) , (3.5) 

Ct^ ot cosh T 

so schematically 

(A(a, T, 0,(t>)^ I dp S;' {a)Ypirn (r, 9, 0) . (3.6) 



Here YpimiT, 9, (f) — fpi{T)Yira{9, (f). The Yim{9, (f) are the usual spherical harmonics. Since all the mass dependence 
is in Sp [a), through 



i/2 



(3.7) 



consider a dependence first. 

The equation of motion for Sp {a) is 



+ 3i/ cot Ha-^ + "^^Y/^ - m\a)]S; (a) = 0. (3.8) 
oa^ da sm Ha 

The most convenient basis, different from the one usually used to describe the tunneling as in is the (u,t) 
coordinate system used to describe the instanton in RecaUing that cos[Ha] — tanhu, 

S;[cos[Ha]) = (3.9) 
^ a(u) 

where 
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Finn) = ^{a.-+r(l - ^p)P+;nUnhu) - ^^a,'-T{l + ^p)P;;'^{tanhu)} 



-x—tan 



h« j^{a.'+(i±f )^2Fi(-z.', i.' + 1; 1 - ip; ^) (3.10) 



lg|a.,_(i±f )^2i^i(-;^', i^' + 1; 1 + i^)} 



The coefficients 



depend on momentum, 



= .i±vEm^ ,,n, 



C1(P) = J-rfl + IS^fi^l 

r (r^ — r(l-zp)r(zp-iy') sin[Ku']sin[T!{v' -ip)] in -, r,\ 

_ 2^2 r(l-zp) sin-n-t/' ]_ 



T{l+ip) sinh'^ 7rp r( — ip— i/')r(l— ip+i/') 

The coefficient Ci(p) is real for values of v' considered here, C2{-~p) — C-zip) finip^oo.i^' fixed C2{p) = 0. 
The functions F^^ are normahzed as 

duF:^iu)Flf'*iu) = ^^Sip-p') . (3.13) 

There is a continuum of modes for < p < oo. (The r dependence introduces some more p dependence so that the 
fuU wavefunctions are nonsingular at p = 0.) 

For the special cases of i^' = (conformal mass, = 2H^) and i^' = 1 {AP = as is the case after tunnchng), 
one has 

rpip _ e'P" 
-'^0 



4^VlPl ("5 1 4^1 

pip _ e'^" jtanhu-ip) \^-^^J 

1 ~ 47r^/H 



When < I/' < 1 (ie when AP/2H'^ < 1), there is an additional normalizable state which was shown to be part 
of the basis of states in They showed that it is necessary to include this state in order to obtain the proper 

Wightman function for the Bunch-Davies vacuum of de Sitter space. 

The Wightman function is the sum over positive frequency modes times the negative frequency modes (a vacuum 
dependent assignment), 

G+{x,x') = {0\^{x)(t){x')\0) 

= (0| Efc,fc'{("+,fc(2;)a(fc) + u-^kix)a^k))} , . 

x{{u+,k'{x')a{k') + u^,k'{x')aHk'))m 

The modes u±^k{x) obey the Klein Gordon equation and are normalized according to 

{u+,k{x),u+,k'{x)) ^5{k-k') 

(u_,fc(x),M_,fc-(a;)) =-5{k^k') (3.16) 
{u+^k{x),u^^k'{x)) =0. 

where the inner product is the Klein Gordon inner product: 

(0i(a;),02(a;))--i / d^^^[-g^{x)]^ {<j>^{x)d^,(^l{x) ~ d^c^^{x)cj,l{x)} . (3.17) 



The integration is over E, a (globally) Cauchy surface, i.e., a spacelike hypersurface which every non-spacelike curve 
intersects exactly once p2| . 

The additional normalizable state for Q < v' <\ has p — iv' , 

F-y ^ — 3_(coshu)-'^' Q<v'<l (3.18) 
47r,/|p| 



6 



It is for a discrete value of p and hence corresponds to a state which normaHzes to a constant times a Kronecker delta 
function, not a Dirac i5 function. 

There is a discussion of the properties of these states in . When continued to the inside of the light cone of the 
bubble center, region I, these states appear as supercurvature modes which are naively unnormalizable. As explained 
in |3^, this lack of normalizability is due to the curves of fixed time T failing to be Cauchy surfaces. (The light cone 
of the bubble center, a non-spacelike curve, is T = 0, and never intersects curves with T ^ 0.) They show that for the 
other modes, with > 0, the inner products on curves of fixed time T coincide with their norms on Cauchy surfaces 
because of sufficiently fast falloff at infinity. The supercurvature modes are normalizable on a Cauchy surface, for 
example along the t = surface outside the forward light cone , and the result is independent of Cauchy surface 
as it should be. 

The i = mode for i^' — > 1 has zero norm, as shown in , perhaps due to overlap with a translational zero mode 
in the system, similar to what happens in ]T^ . Aside from this case, the normalizable discrete state with special 
p = iv' occurs as part of the Bunch-Davies vacuum for all £, m, when < j^' < 1. 

When r dependence is included, the Wightman function can be used, via the Klein Gordon inner product, to project 
upon positive or negative frequency modes for a given set of basis functions. This was used in P,|20|| to identify the 
£ = positive freqency mode functions for the Bunch Davies vacuum 



^+{u,r)=HcoshuF^;P{u)^-^l^Z=^J-^ < p < oo (3.19) 

y'e^P — e~^P cosh r 

with corresponding negative frequency modes 



q>Zpiu:r)=HcoshuF^'Piu)^ < p < oo . (3.20) 

Ve^P — e~^P coshr 

For general £, the r dependent factor fpi (r) is proportional to (from jSl] ) : 

+ l .f^„ J^(-'P) Pn-.Vi'-l' \' sinp(T+i7r/2) 

— H J-J Y '^^^(-ip+i+l) ^ Usinhr/ coshT 

For simplicity, the discussion will be restricted hereon to ^ = m = 0, since the matching across the bubble wall 
depends on a and not r. So, suppressing £, m dependence, the operator describing fluctuations is 

5^{u,t) = dp{(t)+{u,T)a{p) +(ptp{u,T)a{-p) + h.c.} ,^ 22) 

+ Ep^„>p,„, ^)a(p*sc) + /i.e. 



with the normalized in the Klein-Gordon inner product, eqn (3.17) 



IV. INITIAL CONDITIONS, MATCHING ACROSS WALL AND NORMALIZATION 

Given the positive frequency modes in the false vacuum, the next step is to analytically continue to Euclidean time 



and match them on to the initial conditions for fluctuations around the instanton, eq. (2.9). However, the initial 
conditions for the instanton are not directly visible in the (cr, r^;) or {u,te) coordinate system. In these coordinates, 
the bubble wall is present for all time te, so there is no false vacuum at any time te- However, these coordinates are 
useful since a constant value of 6(cr), a{u) deflnes the bubble wall. 

This problem is emphasized in p8[ |, where another coordinate system is suggested which maps (cr, r^;) — > (ct, f^). 
The new time coordinate te coincides with te along the line te = 0, the nucleation time, where analytic continuation 
to Lorentzian signature takes place. The {a, te) coordinate system can be used to specify initial conditions, but it 
is not useful in practice for matching across the bubble wall. At early times te in the (ct, f^) coordinate system, 
the bubble wall has not yet formed, so that the vacuum can be identifled. One can do the matching using either 
coordinate system. The modes corresponding to the false vacuum Bunch-Davies modes in the (tr, t^;) coordinates 
coincide with those in the far past, at early times te, in the coordinate system which has the proper initial conditions 
(false Bunch-Davies vacuum). So, knowing that the (tr, f) coordinate system exists justifies starting with false vacuum 
Bunch-Davies modes exterior to the bubble, as these are continued from modes at early te^ 



I thank Arley Anderson for discussions about this. 
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An important consequence of not specifying the initial conditions at fixed te (th e bubble wall exists for all te, the 



initial condition is without the bubble), that is, using the basis functions eq. ( 3.1E ) for matching, is that the resulting 
wavefunction matched across the bubble wall is not automatically normalized. So the wavefunction is normalized 
after the nuclcation time te — 0. The Minkowski space analogue of this was done in, e.g., |19[|. 



With this understood, the first step is to analytically continue the positive frequency wavefunctions (3.19) of the 
false vacuum to Euclidean time, r —> Ite, 



A,y{u,rE)^c^U^,TE)^HcoshuF^^^{u y . (4.1) 

Ve^P — e~^P cos Te 



Recall ^' = yj — — i, and F^^{u) is given in eq. (3.1C). These are to be matched aX u = uj^ onto linear 
combinations of modes corresponding to the true vacuum, where Af^ w 0, that is v' ^ 1. The basis functions for true 



vacuum modes are 



^ ^^(e— eP^^ -e— e^P^^) 



B±p{u, te) = Ht coshju + 5)F^'P{u + 5)'-^ , , '- . (4.2) 

2Ve^|p| - e-'^IPl COSTS 



where F^{u) is given in eq. (3.14) 



Because the te dependence is continuous, one can only match modes with equal values of |p| (the matching has 
to be true for all times, so the time dependence has to be the same on both sides). Both B±p(u,te) are positive 
frequency modes, so this does not look like particle creation. Continued inside the light cone, where the space and 
time dependence get interchanged, it can be interpreted as causing particle creation. Matching the functions and 
their first derivatives at u = ur, 

Ap{u,TE) = apBp{u,TE) + f3pB^p{u,TE) ,^ g-, 

duAp{u,TE) ^ apduBp{u,TE) + ^pduB^p{u,TE) , 



gives 



evaluated at u — ur. 



o _ -3..-Bp/p+Bp3..Ap (4.4J 
PP ~ B„a^B-„-B-„d^B„ 



So for false vacuum mass M — 2H , (using the explicit form of FqP(u) in eq. (3.14)), 



cosh 5—ip sinh(2'Uj^+(5) i 
2p{p—i) L'^^ ^ cosh cosh(ii/i' -|-(5) J 
_ eip(2^R+^) cosh 5-ip sinh 6 
'"^P 2p{p--\-i) cosh Ur cosh(ii 

For more general masses in the false vacuum, the Bogoliubov coefficients are: 



(4.5) 



^py = §(^{(l+p'-(f^' + 2)tanhMfl(tanh(Mfl + 5)+ip))47r/H^J(^i?J 

+ (i/' -ip+ l)(tanh(ui^ + S)+ ip)a^-+T{l - ip)Plf^^{tanhuR) (4.6) 
-{ly' + ip+ l)(tanh(uii + 6)+ ip) a^'_r(l + ip)P~'P-^{tanhuR)} 



and 



f^P.-' = 1p(p+? ^(^+P' " + 2) tanhuR{tanh{uR + S)- ip))4n^\FlP{uR) 

+ -ip+ l)(tanh(ufl + S)~ ip)a^'+T{l - ip)Pj^i(tanhMfl) (4.7) 
-(i/' + ip+l)(tanh(u;j^ + <5) - ip)-^^a^>^T{l + ip)P^,'_^^{tanhuR)} . 

One can verify that as the radius goes to zero, m ^ oo, (5 ^ 0, and these reduce to the expressions of P,pO|, 

C2 

«p ^R^o 1 Pp ^R^o ~"rp7i^i''- (4-8) 

I02I 

In showing this it is helpful to recall that lim^^oo r(l — ip)Pj^f (tanhw) ^ e'P", so that terms cancel between F*? and 
Pj^f|_-^. In addition, as p ^ 00 these go over to the conformally coupled = 2H^ case. 
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apy^l,Ppy^Q. (4.9) 

That is, at short distances, the dependence upon the scales set by the mass and radius drops out. 

These Bogohubov coefficients are for p real. When < i^' < 1, as mentioned earlier, there are also states corre- 
sponding to supercurvature modes in the false vacuum. These supercurvature modes have r dependence which match 
onto eigenfunctions of the Laplacian with u dependence 

P^' (ni -I- A^l — p-v'{u+S) (tanh(M+a)+iy') 

or F--' {u + 5) = e'''(-+^) ) ^ ^ • 

non-normalizable on any surface. The first eigenfunction is well behaved as m — > cxd and can be connected across the 
wall at M = ur to eigenfunctions which are well behaved at u = — oo, however, it is not possible generically to match 
the original vacuum mode and its first derivative (two constraints) across the wall without the second eigenfunction 
which blows up at u/? = oo. The power spectrum inside the bubble may be computed anyhow, e.g. as is shown in 
for general supercurvature modes, but it is possible that the lack of normalizability means that there is some 
instability (and thus inconsistency in this naive treatment). The way out of this is p4| to not include the state which 
has become unnormalizable due to the wall, but to look for any new states that may have become normalizable at the 
same time, again due to matching across the wall. Thus one has an added degree of freedom (the value of p) which 
can be used to try to satisfy the matching constraints. As they show, a normalizable supercurvature state is present 
for M^/7J^ < 1 for i? = 0, and for larger values of M when i? ^ 0. 

Similarly, there is also a supercurvature mode present in the 'true' vacuum of the bubble interior, where v' \. It 
does not have the same time dependence as any state in the false vacuum, and so is not excited by the continuation 
of the false vacuum fluctuations across the wall. It is present in some cases, e.g. if there is no bubble, |^,^, or if 
the field providing density fluctuations is massless on both sides of the bubble wall |^ . If the mass of the fleld in the 
false vacuum is large enough, the question of these supercurvature modes does not occur. These modes will not be 
considered further here. 

Once the modes have been matched across the wall, the wavefunctions are analytically continued to Lorentzian 
signature at = = —it to give the wavefunctions for a bubble outside of the light cone of the center of the bubble. 
This is straightforward to do by substituting t = ite- 

The resulting wavefunctions are 

'?!'+p,unnorm("7T) = 9{ur - u)Apy {u,t) + e{u - UB){apBp{u,T) + l3pB^p{u,T)) , (4.11) 

where 0{u) is a step function. As mentioned before, (j)~^p unnorm(''^: ''') unnormalized. In fact, it is not even|^ orthogonal 
to 0!tp.unnoim("j ''')■ This cau be remedied by using the same method in [po[ . 

As the T dependence is the same for both fields, and corresponds to a properly normalized wavefunction (as would 
be automatic if there were no bubble wall, which also has the same r dependence), it can be separated out first. 



a{u) 2Ve'^\p\ - e-'^IPl cosr^ 



+p,unnorm 



^u,r) = ~ : — ^ . (4.12) 



Then, the T'^p{u) are eigenfunctions of a self-adjoint operator with eigenvalue depending on p^, and so eigenfunc- 
tions with different eigenvalues are orthogonal, 

duT^P(u)T^P (u) = ^'5(p + p') + ^^'5(p - p') (4.13) 

The coefficients Di{p), D2{p) can be read off from the asymptotic behavior at u — s- ±00. After some algebra they are 
seen to be: 

Diip) = HiapP + + 1] = A(-P) 



D2{P) 

Then the normalized wavefunction is 



C2 

Pf-^P ' 2Ci 



I thank M. Bucher for pointing this out and for suggesting the method used below. 
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(4.15) 



where 



b±,p{uR) 



\D, 



i±^i-\D2\yDi 



(4.16) 



This normaUzation can be checked by a different calculation once the modes are continued into the future light cone 
of the bubble center, as described below. The other method of normalization has no integration over (the analytic 
continuation of) u, and so does not depend on continuity the same way. The two calculations agree. So the operator 
describing fluctuations, 



is as in eq. ( p.22| ), with (j)p{u, r) given in eq. ( [4.15[ ) 



V. CONTINUATION INTO OPEN UNIVERSE, POWER SPECTRUM 

The operator for fluctuations is next continued inside region I to find its effect in the open universe. For u > ur, 
approaching the light cone, 



0p(w,r)|„>„„ = {b+.p{uR)ap- ^j^b^^p{uR)P^p)Bp{u,T) 
+ {b+^p{uR)Pp - ^jj^b^^p{uR)a^p)B^p{u,T) 



(5.1) 



At the light cone, cr — )■ (u — > oo). Inside region 1, a — —iT, t — x — i7r/2, and time and space get interchanged. The 
analogue of u in region II is the conformal coordinate ry, obeying = tanh(T/2). 

The a, t coordinate system becomes singular at the light cone. The correspondence of wavefunctions across the 
light cone was found in |^ , and can be taken over directly, as they emphasize, since only the short distance behavior 
of the fields is used: 



e-*P"coshu^4^^ -> 

cosn r 

e*P" coshu-£^ ^ 

cosn r 



(5.2) 



e^P^COShu^^ ^ (_i)4H£Xg-(*p+l)r, 
cosn T ^ ' Sinn v 



Only half of the modes propagate inside the light cone for massless fields, the other propagate along the light cone 
instead 1^. Since the matching depends only on the coefficients of ir in the exponential, these rules can be used 
directly to match the modes for general i if desired. 

The basis functions then become: (one can think of u here as u + S in the earlier basis functions) 



Bp{u,r) ^„ -,i7^^^=^=e(-P-i)'™ ^ B,(^,x) 



R (y r) -iHrr— e 2 (ip-Dy sinpx ^ n / \ 



(5.3) 



With this continuation across the future light cone of the bubble center, eq. (3.22) becomes 



H{V,X) = dp {b+,p{uR)ap ~ ^j^b-^p{uR)f3^p)Bp{r],x) 

+ {b+,p{uR)(3p- ^^b^^p{uR)a-p)i3^p{ri,x) a+{+p) 

+ [{b+,p{uR)a^p - £ft^b.,p{uR)Pp)B.p{rj, x) ^^'"^^ 

+ {h+,p{uR)(3-p- ^f^b-^p{uR)ap)i3p{'n,x) a+{-~p) 
+h.c. 

(Recall there are no discrete states for the values of AP being considered.) 

This analytically continued wavefunction can be verified to have proper normalization, providing a check on the 
calculation of normalization given earlier. The future in de Sitter space has two light cones which must be included, 
region I, where the open universe interior to the bubble is found, and that of the antipode of the center of the bubble 
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(for descriptions see, e.g. P, |3l] , |3^ ), region III. (These are regions L,R in the notation of |3l[].) For our purposes, 
continuing into region III can be done by using symmetry and then verifying that a normahzed wavefunction in 
region II remains normahzed when continued into region /, /// using these rules. The modes continued into region 
III are exterior to the bubble wall, and consequently have mass corresponding to the false vacuum mass. As time 
and space are interchanged in these forward light cones, the Klein Gordon inner product now involves Wronskians 
of the analytically continued u dependent part of the wavefunctions. The inner product is calculated on a fixed 
time surface in both forward light cones and then added. Although, as stated before, the fixed time surfaces are not 
Cauchy surfaces, the wavefunctions under consideration are the same ones (restricted to each light cone) discussed in 
pl[ . They do not correspond to supercurvature modes and consequently have sufficiently fast falloff at infinity that 
their inner products on the fixed time surface inside the light cones correspond to their inner products on a Cauchy 
surface. It is also usef ul t o know, in order to show this, that — \f3p\'^ — a^,^ — a^,_ which can be shown with 



the definitions in eq. (4.4) and properties of Wronskians. 

The fluctuation S(f) inside the light cone sources gravity, and can be matched directly onto fluctuations of the 
gravitational potential <& using results of as follows. Since the fluctuations are small, they can be described by 
expanding the metric around de Sitter space and the field (j) around its classical value, as was done in also see 
p2| for another discussion. The result is that a source (j) « e^^'^P~^^^ corresponds to a "gauge invariant gravitational 
potential" 

$ « l^L^e^i^f+i)" . (5.5) 
±ip+2 ^ ' 

Then, as |^ show, when H remains constant and the potential is linear, this matches onto the exact solution: 

$ = ^^,i±^P+i}n[^ _ P±l-e^^] . (5.6) 
±ip + 2 L 3(pT«) 

For the 6(j) found here, the corresponding positive frequency $ is: 

^(n — j -^'^'^^-*' rV\ f°° dp sinpx 1 

+e^{b+,,{un)f3p - £0b^^,{un)a^p)^e^^^l - ^e-^'])a{p) (5.7) 

+e¥(6+,,(u«)/3_, - £ft^b.A^n)ap)^e~^PV[i „ _£_i_e2.])a(-p)} 
The power is defined by taking 

($(x = 0, 77)$(x, ry)) = / dpp — - — P<;>(p,ri) + discrete states (5-8) 
Jo smhx 

and isolating the growing mode by taking rj 0^{T oo). This gives 

p(l+p^)2sinh7rp 9H'i, 

{\e^{b+AuH)ap- ^^b-^URjfi-p) 
-^^,e^{b+Aun)Pp - ^^b.,,{un)a^,r 

+ \{e^{b+,j,{uR)a^p - ^^b^AuR)P,) (5.9) 
-H^e¥(6+,,(u«)/3_, - £2t^b^,AuR)a,)\'' 

= {GV^^Y p(l+p^)2 sinh Trp [^''^ + ^^'^^ 



Here, the identity 



= i(l + \a,\' + |/5pp - f ap/3p - a_p/3_p - |g^p) 



was used. 
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Taking i? — s- 0, i.e. ap 



this reduces to the case found in l2C 



9iJ2 p(p2 _^ i)2sinh7rp 



+ e 



C2 p + i C2 p — i 



Cip- 



(5.11) 



Generally < 1, so in the i? = case, this means that the p dependence is in an envelope between coth ^ and 



tanh ^ , going to {p(p' 



1))^^, the scale invariant spectrum, approximately when p>2. 
For < p < 2 a wide range of values are possible, depending upon the specific values of 5, R, M/H. For p > 2, 
when the radius is nonzero, the last term (the last line) can only be significant if it grows faster than the prefactor 
sinhTrp"^ goes to zero, in such a way as to be comparable to the first term e^^ > 250. 
The magnitude of this term can be bounded as follows. It can be rewritten as 



4i?e[e 



i4,P + i 



yA- 



P - I (71^(2/2 - 1) + 2/2 + i)A - ycM' + 1) 



(5.12) 



where 



A 



y 



IC2I |/3_p| 



(5.13) 



and it has been used that japp 

_ /3p a-p 



- = V^~^c and that = ^^ec- (Also recall that a_p = a*p = /?*.) The 
phase d is real, e*'^ = ; and will not be of interest here since the magnitude will be shown to be bounded and 

small. From their definitions above, one has 0<e<l, 1 < y and |A| = 1. 

It is possible to extremize the above expression in A and y because the extrema of a ratio of two quadratic 
polynomials is given by a quadratic equation and can be solved analytically. However, this is of limited use. In the 
space spanned by physical values of ec, y, A, the expression was not found to ever be larger than two in absolute value 
(fixing one parameter, plotting for the full range of the other two, and changing the parameter that was originially 
fixed). In the limiting value of Ec = 1, the magnitude can be seen directly to be two and when — s- it is of magnitude 
less than two. For any mass value, as p is changed, one runs through a set of allowed values of Cc, y, A, but for p > 2 
this term is effectively irrelevant, and the power spectrum becomes scale invariant to a high degree of accuracy. That 
is, if equation 5.12 is less than or equal to two in magnitude, then the quantity in brackets in equation 5.9 is between 



coth7rp/2 and tanh7rp/2, just as in the zero radius case. 

The effects from the vacuum fluctuations propagating through the bubble wall must be combined with other 
fluctuations caused by the presence of the wall. In quantizing the field theory around the classical bubble solution 
4>b{'^) in flat space, there are two special modes arising, which remain even in the thin wall limit [ pT| . One is a 
translational zero mode, d^4'b{(T), where = dx^, and and one is a mode coming from sh iftin g the radius, da(t>bicr). 



The latter state was identified by |ll|] as having p^ = —4 in the basis functions of eq. (3_^), and appears also in 
the de Sitter case if one assumes that H changes negligibly during tunneling The former corresponds to 

(nonorthogonalized) ^ = 0, 1 modes for this value of p. As the explicit form of both is related to derivatives of the 
classical background, which only varies significantly on or near the wall, they tend to be localized on the wall, and 
their degrees of freedom can be related to fluctuations of the wall. Although localized, fluctuations in the wall can, 
at later times, change the definition of where the 'time' is zero. So these fiuctuations can produce an asymmetry in 
space now, since different regions would be at slightly different times in their evolutions. In 0,^,^ 12 1 the bubble wall 



fiuctuations are discussed and contributions to asymmetries and to temperature anisotropics in the CMB calculated. 



VI. CONCLUSIONS 

In this paper the results of [^,0 were extended to include a finite wall radius and small vacuum energy difference 
before and after tunneling. The power spectrum appears to generally go over at p ~ 2 to an approximately scale 
invariant spectrum just as in the zero radius case. As most of the discriminating effects found here are at very large 
scales, cosmic variance may limit their observability. The scale invariant spectrum that seems to arise at smaller 
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scales (no exceptions were found and perhaps it can be shown rigorouslyQ) means that varying the parameters R, S 
and false vacuum mass M over a large range for thin wall bubble models gives a generic prediction. 

For a given model, these effects can be translated into effects in the cosmic microwave background by combining 
them with the fluctuations from the wall modes and a model for matter such as cold dark matter. Modifications of 
this simplified model (wall thickness, tilt in the final spectrum) should be considered, again to quantify the sensitivity 
of the models to variations in parameters. Also, so far, the gravity wave perturbations have only been calculated for 
the wall fluctuations themselves ||l2| . 

Along other lines, it would be interesting to see what sorts of potentials naturally arise in particle theory models. 
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